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Abstract — Recent work by Divsalar et al. has shown that 
properly designed protograph-based low-density parity-check 
(LDPC) codes may have minimum distance linearly increasing 
with block length. This fact rests on ensemble arguments over 
all possible expansions of the base protograph. However, when 
implementation complexity is considered, the expansions are 
frequently selected from a smaller class of orderly expansions. 
For example, protograph expansion by cyclically shifting con- 
nections generates a quasi-cyclic (QC) code. Other recent work 
by Smarandache and Vontobel has provided upper bounds on 
the minimum distance of QC codes. In this paper, these bounds 
are expanded upon to cover puncturing and tightened in several 
specific cases. We then evaluate these upper bounds for the family 
of protograph codes known as AR4JA that have been proposed 
for deep space usage in the CCSDS experimental standard. 
Finally, we note that at block lengths more than 4400 bits these 
upper bounds fall below the ensemble statistical lower bounds of 
Divsalar, generating a new perspective on the minimum distance 
for the AR4JA codes in the standard. 

Index Terms — binary codes, block codes, error correction 
codes, linear code, sparse matrices 

I. Introduction 

Avery important class of modern codes, the low-density 
parity-check (LDPC) codes, had their start in the seminal 
work by R. Gallager fTI about 50 years ago. Work on these 
codes remained mostly dormant for decades, with the impor- 
tant exception of R. M. Tanner's work on graph-based code 
construction Q. Properly designed LDPC codes exhibit very 
low SNR thresholds in their error rate performance. However, 
there has been a tradeoff evident between SNR threshold and 
error floor performance. One technique to lower error floors 
in LDPC codes reduces the number of small cycles in the 
graph and optimizes the local multiplicities of the small cycles 
||3). Similarly, the ACE algorithm [4] for placing edges in a 
graph-based code brings down the error floor substantially by 
preventing small cycles from clustering on low-degree variable 
nodes. 

Another important property limiting the error floor of any 
code is the minimum distance, yet relatively little work on 
designing LDPC codes for large minimum distance has been 
reported. The minimum distance is also important in under- 
standing the likelihood of undetectable error patterns, a critical 
parameter in certain applications such as data storage. 
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Fig. 1. Simple protograph and corresponding protomatrix A. 

Most codes within the ensembles of codes based on pro- 
tographs with certain properties have been shown to have 
minimum distance linearly increasing with block length (5), 
(6). These protographs together with the ACE algorithm have 
been used to design LDPC codes for deep space usage in the 
CCSDS experimental standard j7). These codes, as specified, 
are punctured LDPC codes and fall into the class of Quasi- 
Cyclic (QC) codes. A separate body of work on QC-LDPC 
codes exists (8), |9), including recent work concerning upper 
bounds on the minimum distance flO) . We attempt to bring 
these works together by extending the bounds to punctured 
QC-LDPC codes and tightening the bounds where possible. 
The authors previously presented a summary | [TT[ of this 
subject, which is treated here more completely. 

II. Protographs and AR4JA 

Protographs were introduced as a way to impart structure to 
the interconnectivity of graph-based codes fl2") . Protographs 
themselves are a subset of the multi-edge type graphs JT3}. 

A protograph is a Tanner graph with a relatively small 
number of nodes, except that parallel edges are permitted. A 
protograph, G = (V, C, E), consists of a set of variable nodes 
V, a set of check nodes C, and a collection of edges E. Each 
edge, e E E, connects a variable node, v e £ V, to a check 
node, c e e C. A useful refinement is to allow the variable 
node set V to contain punctured (i.e., not transmitted) nodes. 

A simple protograph is shown in Fig. [T] with three variable 
nodes, two check nodes, and five edges. The accompanying 
protomatrix fully describes the graph. The entry in the ith 
row and jth column of protomatrix A indicates the number 
of edges connecting the zth check node to the jth variable 
node within the corresponding protograph. The labeling of the 
protograph in Fig. [T] indicates node types. 

The derived graph is constructed by replicating the proto- 
graph many times and interconnecting the copies. The pro- 
tograph code is defined by the resulting derived graph. All 
copies of check node A, are termed "type A" check nodes. 
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Fig. 2. Copying the protograph N = 3 times. 
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Fig. 3. Finally, interconnecting the copies generates the derived graph. 
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Fig. 4. AR4JA protograph, rate-1/2. The transmitted variable nodes are 
shown as solid circles, the punctured variable nodes as outlined circles. 
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Fig. 5. AR4JA protograph, rate-2/3. 
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Similarly, all copies of variable node 1 are termed "type 
1" variable nodes. The interconnection process proceeds by 
treating all copies of an edge as an edge set, and swapping con- 
nections only within each edge set. This rule prevents nodes 
from changing degree and maintains the graph connectivity by 
node type. 

The main advantages of protographs are that degree-one 
variable nodes and punctured variable nodes may be in- 
troduced in a structured way. The optimization of standard 
irregular LDPC codes by density evolution does not allow for 
degree-one variable nodes, but produces a significant fraction 
of degree-two variable nodes. Thus, it is natural to conjecture 
that degree-one variable nodes may bring benefits if they can 
be incorporated into the code fT3) . An additional advantage 
of protographs is that decoder implementation may be less 
complex due to the structured interconnections. 

Figs. [2] and [3] illustrate the process of making N — 3 copies 
of the protograph of Fig. [T] and interconnecting them to gener- 
ate the derived graph. The parity-check matrix corresponding 
to the derived graph of Fig. [3] is shown below, divided into 
submatrices so the relationship to the protomatrix A of Fig. [T] 
is evident. 



protomatrices are shown in Q and Q, respectively. 
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The protograph for the rate r = 1/2 AR4JA code |31 
is shown in Fig. [4] We follow the convention of showing 
transmitted variable nodes as solid circles and the punctured 
variable nodes as outlined circles. The protograph of the rate- 
1/2 code is extended to rate-2/3 by adding two degree- 
four variable nodes as shown in Fig. [5] The corresponding 
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The variable nodes have been numbered in the figures to 
correspond to columns of the protomatrices, from left to right. 

The AR4JA family of protographs continues to increase the 
offered code rate options by adding more pairs of degree-4 
variable nodes. As the protomatrix Q shows, the rate-4/5 
protograph has 11 variable nodes altogether. In all cases, the 
variable nodes corresponding to the right-most column of the 
protomatrix are punctured, i.e., the variable nodes of degree-6. 



A r= 4/ 5 — 



00000000102 
3 13 13 1110 13 
1 3 1 3 1 3 1 2 2 1 



(4) 



The name AR4JA is derived from the operations indicated in 
Figs. |4] and [5] and signifies a resemblance to the Accumulate- 
Repeat-Accumulate (ARA) codes in protograph form. For 
AR4JA a partial precoding by accumulation ("A") is followed 
by a repetition by 4 ("R4") and its design ends with a jagged 
accumulation ("JA"). While a standard accumulation stage 
contains solely degree-two variable nodes, this accumulation 
stage, with its extra edge, appears "jagged." 

Divsalar et al. presented techniques for calculating the 
asymptotic ensemble weight enumerators for protograph-based 
codes 0, (6j, |14| . From the derived expression of the weight 
spectrum, the typical minimum distance ratio <S mm can be 
found, if it exists. Then, with high probability the minimum 
distance of most codes in the ensemble increases linearly with 
block length n with proportionality constant (5 m ; n . The AR4JA 
rate-1/2 protomatrix (|2| has 5 m \ n — 0.015. The switch from a 
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standard accumulation stage to the jagged accumulation stage, 
with its reduced number of degree-2 variable nodes, allows 
the AR4JA protograph to have linearly increasing minimum 
distance (5). 

III. QC Expansion and Polynomial Representation 

The codewords of a block code may be divided into non- 
overlapping subblocks of N consecutive symbols. A quasi- 
cyclic (QC) code is a linear block code having the property 
that applying identical circular shifts to every subblock of a 
codeword yields a codeword. QC codes are a generalization 
of cyclic block codes. Thus, a conventional cyclic block code 
is QC with a single subblock. 

A binary QC-LDPC code of length n = LN can be 
described by an m x n sparse parity-check matrix, H £ F™ xn , 
with m = JN. The code can also be described in polyno- 
mial form, since there exists an isomorphism between the 
ring of N x N circulant binary matrices and the ring of 
polynomials with binary coefficients of degree less than N, 
F2[x]/ (x N — l). Addition and multiplication in the latter ring 
follow standard polynomial operations with coefficients in F2, 
modulo x N — 1. (All the rings in this work are commutative 
rings containing a multiplicative identity.) 

A right circulant matrix is a square matrix with each 
successive row right-shifted circularly one position relative to 
the row above. Hence, circulant matrices can be completely 
described by a single row or column. We use the left-most 
column convention to define the isomorphism as in |10|. 
Hence, the isomorphism maps the set of N x N binary 
circulant matrices to polynomials by assigning to polynomial 
coefficients of increasing order the left-most column entries 
of the circulant matrix from top to bottom. The polynomial 1 
maps under the isomorphism to the N x N identity matrix. 
A few examples of the mapping (indicated by H>) for N = 3 
are shown below. 
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This isomorphism requires that care be taken when repre- 
senting multiplication of a circulant matrix M by a binary 
row-vector v = {vq, V\, . . . U/v-i)- We can associate the 
polynomial M(x) with the matrix using the technique just 
described, and associate v(x) = vq + v\x + . . .vn-ix 1 ^^ 1 
with the vector v. The product (Mv T ) = vM T maps to 
the polynomial M(x)v(x) modulo x N — 1, and the product 
vM maps to the polynomial M (x~ 1 )v(x) modulo x — 1. 

Our application of this isomorphism will require that we 
define the weight of elements in the polynomial ring. Given 
a(x) £ ¥2{x]/ (x N — l), we define its weight wt(a(x)) to 
be the number of nonzero coefficients in a(x). Thus, the 
weight of the polynomial wt(a(x)) is equal to the Hamming 
weight of the corresponding binary vector wjj(a). Addition- 
ally, for a length- N vector of elements in the ring, a(x) = 
(ao(x),ai(x), . . . ,cln-i(x)), we define its weight to be the 
sum of its component weights, wt(a(x)) = X^o* wt(dj(x)). 



As we will be using elements from a ring and not a field, 
there are a few important concepts to keep in mind. The 
elements in a ring are not necessarily invertible, and the ones 
that are invertible are called units. Also, a ring may include 
zero divisors. A zero divisor or factor of zero in the ring R is 
a nonzero element a £ R, such that ab = for some b £ R, 
b 0. For a simple example, consider R = Z/6Z, the ring of 
integers modulo 6, with zero divisors a = 2 and 6 = 3. Units 
in the ring cannot be zero divisors. 

The elements of weight 1 in the polynomial ring 
F2[x]/ (x N — l) are the monomials, which are units in the 
ring. Specifically, the inverse of the monomial a(x) = x % , 
< i < N, is the monomial [a(x)] _1 = x N ~ l . The 
monomials in the polynomial ring map to cyclic permutation 
matrices under the isomorphism. A cyclic permutation matrix 
is a permutation matrix (a square matrix of ones and zeros, 
such that each row and column contains a single one) that is 
also circulant. 

The nonzero elements with even weight in the polynomial 
ring Fa[x]/ (x ~ l) are zero divisors. For instance, let o(x) 



be the all ones polynomial x 



N-l 



„N-2 



1. Then, 



for any even weight polynomial b(x) in the ring, b(x)o(x) 
modulo x N — 1 is zero, for any positive integer N. Odd weight 
polynomials may sometimes be zero divisors, such as x 3 +x+l 
in the ring F 2 [x]/ (a; 7 — l). 

As we are interested in the connection between protographs 
and QC-LDPC codes, we focus on parity-check matrices, H, 
that are in J x L block matrix form, in which 

Ho.o • • • H ,l-i 

H = ; ••. ; 

Hj-1,0 • • • Hj_i,£_i 

where each submatrix H,-^ is an N x N binary circulant 
matrix. Let hj t i t i t k £ F2 be the entry in the lih row and fcth 
column of H^j. Then, using the left-most column convention, 
we may state EL,-.; = J^^q 1 hj&sfif** where I s is the N x N 
identity matrix circularly left-shifted by s positions. Now, 
we can write the polynomial parity-check matrix, H(x) £ 



[F 2 [x]/(^-l>] 



JxL 



H(x) 



hofi(x) 



h 



j-i 



o(x) 



h ,L-i(x) 
hj-ix-i(x) 



where h hl (x) = J2 s =a fy?,i,s,oX s . 

Furthermore, we will be interested in the weight of each 
polynomial entry of H(x) (or equivalently, the row sum of 
each submatrix of H). The JxL weight matrix of nonnega- 
tive integers is 



wt (H(sb)) = 



wt (ho i0 (x)) 



wt (/ij_i i0 (x)) 



wt (h ,L-i(x)) 



wt (hj-i t L-i(x)) 



Note that the protograph connection is evident as the QC- 
LDPC weight matrix wt (H(x)) is exactly the protomatrix, 
A = wt(H(x)). 
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Just as the matrices used to describe QC codes are conve- 
nient in polynomial form, so are the codeword vectors. Define 
a vector of polynomials c(x) G [F 2 [a;]/ (x N — l)] to be 

c(x) = (co(x),ct(x), . . .,c L -i{x)) . 

The elements Cj(x) of the codeword polynomial vector c(x) 
correspond to the codeword subblocks defined at the start of 
this section. To define the set of codewords for a QC code, we 
note that the set of all c such that Hc T = T (with elements 
in F2) maps to the set of all c(x) such that H(x)c(x) T = T 
(with elements in the ring Fajo:]/ (x N — l)). 



IV. Minimum Distance Bounds for QC Codes 

In this section we review the Hamming distance upper 
bounds of [10 1 and extend them to punctured versions of QC- 
LDPC codes. 

We will use the shorthand notation [L] to indicate the set of 
L consecutive integers, {0,1,2, ... ,L — 1}. We will use the 
common backslash notation to exclude a member from a set. 
For example, the set S\i contains all the elements of S except 
element i. Additionally, A5 is the submatrix of A containing 
the columns indicated by the set S. Similarly, the vector 5 
is the subvector of vector a 3 containing the elements indicated 
by the set S. By default, we use row vectors throughout. 

The permanent of the J x J matrix B = [bji] is defined to 

be 



perm(B) 4 £ JJ 



(5) 



where the summation is over all J! permutations a of the 
set [J], and <j(j) is the jth entry of the permuted set. The 
permanent is similar to the determinant of a square matrix, 
which is 

det(B)4^ s ign((7) J] b jMj) , (6) 
ie[J1 

where sign(cr) equals +1 if a is an even permutation and —1 
if a is an odd permutation. When the elements of B belong to 
a ring of characteristic two, perm(B) = dct(B), as addition 
and subtraction are interchangeable. Like the determinant, the 
permanent may be computed recursively by taking the cofactor 
expansion along any row or column. For example, the cofactor 
expansion of the permanent of matrix B along the jth row is 



perm(B) 



E 



bj i ■ perm(B( 



[J]\V 



(7) 



where B' is the submatrix of B with the jth row removed, 
for any j G [J]. 



Lemma 1 (Lemma 6 fit))). Let C be a QC code with poly- 
nomial parity-check matrix H(.x) G [Fa [a;]/ (x —1)1 
Let S be an arbitrary size-{J + 1) subset of [L] and let 
c(x) € [F 2 [a;]/ (x — l)] be a length-L vector whose 
elements are given by 



Ci(x) = 



perm (H S \i(x)) ifieS 







otherwise. 



(8) 



Then c(x) is a codeword in C. 



Proof: Let the jth row of H(x) be hj(x). Then, 

hj(x)c(x) T = ^2 h 3 ,i(x) ■ Ci(x) 
ie[L] 

= E h iA x ) P erm ( H 5\i(^)) 



perm 
det 



h. j>s (x) 
H 5 (x) 



(9) 



QVjG[J]. 



The second line is the cofactor expansion of the third line. 
As the elements belong to a ring of characteristic two, the 
determinant equals the permanent. The determinant shown 
must be zero as it contains a repeated row. Since every row 
of H(x) has zero inner product with c(x), H(x)c(x) T = T 
and c(x) is a codeword in C. ■ 
We will require a variant of the min function, which 
normally returns the minimum value of a collection of values. 
For a finite collection of nonnegative real numbers 1Z, let 
1Z + C 72. be the subset of positive elements of 1Z. We define 



mm At = 



min7e+ if TZ+ ^ 
+00 if K+ = 



(10) 



We now begin to bound the minimum Hamming distance 
of the QC code C, which is 



AC) = rnin*wt(c(ai)). 
c(x)ec 



(11) 



Our focus is on two possible upper bounds to the minimum 
distance. The first upper bound will use the code's polynomial 
parity-check matrix to produce low-weight codewords using 
Lemma [T] We will generate as many codewords as possible 
and apply ( [TT) to achieve an upper bound on the minimum 
distance. 

Theorem 1 (Theorem 7 [10]). Let C be a QC code with poly- 
nomial parity-check matrix H(x) G [F2[ie]/ (x n — 1)1 
Then the minimum Hamming distance of C satisfies the upper 
bound 

(C) < min* 53wt(perm(H5\i(x))). (12) 



5C[i 
|S| = J + 1 



Proof: Let S be a subset of [L] of size-( J+ 1) and apply 
Lemma [T] to construct a codeword in C. The weight of the 
resulting codeword c(x) is 



wt (c(x)) 



E 

ie[L] 

E 



wt (Ci(x)) 



wt (perm (H s \j(x))). 



(13) 



Combining (Hi and (13 1, yields (12i as an upper bound on 



the minimum distance as only a subset of codewords can 
be generated by Lemma [T] For some choices of set S, the 
construction in Lemma [T] will yield the all-zero codeword, so 
we must exclude those sets from the calculation of the upper 



bound (12 1. The min* operator does this. 
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Lemma 2. Let B(x) be a J x J matrix with elements in 
the polynomial ring F 2 [x]/ (x N — l). Then the weight of the 
permanent of B satisfies the upper bound 

wt (perm(B(x))) < perm (wt(B(x))) . (14) 

Proof: Let polynomials a(x) and b(x) be in the ring 
F 2 [x]/(x JV - l). We know that wt[a(x)+6(x)] < wt(a(x)) + 
wt(b(x)) as the maximum number of nonzero coefficients in 
the sum a{x) + b(x) is wt(a(x)) +wt(b(x)). Indeed, addition 
in this ring may eliminate pairs of terms. 

Similarly, we know that wt[a(x) ■ b(x)\ < wt(a(x)) • 
wt(6(x)) as the maximum number of nonzero coefficients 
in the product a(x) ■ b(x) is wt(a(x)) ■ wt(6(x)). Note that 
intermediate product terms may be eliminated in pairs as they 
are summed. These properties of the weight permit us to 
complete the proof, since 



parity-check matrix H(x) G [F2[x]/ (s 



wt (perm(B(x))) = wt 



< 



E 



wt 



(15) 



n 6 j>cj) 

^EIl wt ih«u)) 

= perm (wt(B(x))) . 



The following theorem uses Theorem [T] and Lemma [2] to 
create an upper bound on the minimum distance based on the 
code's weight matrix. 



Theorem 2 (Theorem 8 1 10 1). Let C be a QC code with poly- 



JxL 



nomial parity-check matrix H(x) G [F 2 [x]/ (x N — 1^ 
and let A = wt(H(x)). Then the minimum Hamming distance 
of C satisfies the upper bound 

dmin(C) < nun* y^perm (A 5 \ a (x)). (16) 

\s]=j+i ies 

Proof: In part, we combine the results of Theorem [T] 
and Lemma [2] See [10] or Theorem [4] for the proof that the 
potential min* complications have no effect on this bound. 
This potential for complication arises because the S sets 
excluded from (12i by min* will not necessarily be excluded 
from ( fl6l >. ■ 
We now introduce puncturing. By puncturing whole sub- 
blocks of the polynomial codeword vector c(x) we maintain 
the QC property. The subblocks indexed by V C [L] are a 
subset of the L subblocks of c(x). Indices of V may also be 
associated with columns of the Jx L polynomial parity-check 
matrix H(x). 

We begin with an unpunctured QC code C, based upon 
H(x). Next, we define a new QC code, C, by puncturing 
the subblocks of c(x) as indexed by set V . We mark the 
subblocks to be punctured with the symbol "(p" as re-indexing 
would introduce unnecessary notational complexity, and define 
wt(ip) — 0, since punctured symbols are not transmitted. 

Lemma 3. Let C be a punctured QC code constructed by 
puncturing subblocks of QC code C, defined by the polynomial 



: N -l)] JxL . Let the 
subblocks of code C indexed by set V,Vd [L], be punctured. 
Let S be an arbitrary size-{ J+ 1) subset of [L\. Let the length- 
L vector, c'(x) = (cg(x), c^(x), . . . , c' L _ 1 (x)), with c'(x) G 
|F 2 [x]/ (x — l\ , (p\ be defined by 



<{x) 



perm (H 5V (x)) 




ifieS\v 
ifieV 

otherwise. 



(17) 



Then c'(x) is a codeword of the punctured code C. 

Proof: This follows by noting that c'(x) is obtained by 
puncturing subblocks indexed by V from the codeword c(x) 
of Lemma Q] ■ 

Theorem 3. Let C be a punctured QC code constructed 
by puncturing subblocks of QC code C with polynomial 
parity-check matrix H(x) G [F 2 [x]/(x Ar — l)] . Let the 
subblocks of code C indexed by set V, V C [L], be punctured. 
Then the minimum Hamming distance of C satisfies the upper 
bound 



(C) < mm 



SC[L] ' 

si=./+i ies\v 



wt(perm(H 5Vi (x))). (18) 



Proof: Let S be a subset of [L] of size- (J +1), and 
apply Lemma [3] to construct a codeword in C. The weight of 
the resulting codeword is 

wt(c'(x))= ^wt(cUx)) 

= 2^ wt (pcrm(H 5 \j(x))), 

ies\v 

where we use wt(ip) = 0. This expression is combined with 
( [TT] ) to yield our upper bound ( [T8| on the minimum distance 
as only a subset of codewords can be generated by Lemma [3] 

■ 

Care must be taken throughout this work that puncturing 
does not reduce the dimensionality of the code. What we 
mean by dimensionality is the logarithm of the number of 
distinct codewords, which is the rank of the binary generator 
matrix G G F 2 xLJV . Clearly, if a nonzero codeword of C can 
be punctured to the all-zero codeword of C, we have lost 
dimensionality with respect to the original code. 

Lemma 4. Let C be a punctured QC code constructed by 
puncturing subblocks of QC code C, while maintaining the 
dimensionality of code C. Let length- L vector c(x) be a 
codeword of C and c'(x) be a codeword of C obtained by 
puncturing c(x). Then, c^(x) G {0,93} Vi G [L] if and only if 
c(x) = 0. 

Proof: The forward direction follows from the require- 
ment that the dimensionality of the original code be main- 
tained. The reverse direction follows directly from the fact that 
puncturing the all-zero codeword of C produces the all-zero 
codeword of C. ■ 
Since each subblock's contribution to the weight of the 
codeword is nonnegative, the weight of any particular punc- 



tured codeword ( 19 1 must be less than or equal to its weight 
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before puncturing (13) . Moreover, Lemma [4] lets us state that 
the punctured upper bound of Theorem [3] will always be less 
than or equal to the upper bound of Theorem [T] where there 
is no puncturing. 

Theorem 4. Let C be a punctured QC code constructed by 
puncturing subblocks of QC code C, defined by the polynomial 
parity-check matrix H(x) £ [F 2 [x]/ (x — l)] and let 
A = wt (H(x)). Let the subblocks of code C indexed by set V, 
V C [L], be punctured, while maintaining the dimensionality 
of the code. Then the minimum Hamming distance of C 
satisfies the upper bound 



>( c ') ^ E i ,<,nn 

\s\=j+i ies\v 



(A 5V )- 



(20) 



Proof: Let S be a subset of [L] of size-( J+ 1), and apply 
Lemma[3]to construct a codeword in code C, c'(x). From ( 19 1 
we obtain 



wt(c'(x)) 



< 



^2 wt (P erm ( H 5\ l (a;))) 

ies\v 

^2 perm(wt (H^O))) 

ies\v 



(21) 



E 

ieS\V 



perm 



(A 5V ), 



where we invoked Lemma [2] in the second step. 

The next part of the proof justifies the use of the min* 
operator. The potential complication is that for specific choices 
of the set S, namely those which yield the all-zero codeword 
in Lemma [5] the min* operator will exclude their contribution 
to (18), but not necessarily (20 1. In other words, the RHS of 



d2Tb may be a positive value when the LHS is zero. Thus, 
the remainder of this proof assumes c'^x) £ {0, ip} Vz £ [L] 
for a specific choice of S C [L] . By Lemma |4] we know 
that c(x) — for this S. By Lemma [T] we know that every 
J x J submatrix of Hsfi) must have a zero permanent and 
determinant. 

Case 1: If 2~2i£S\v P erin {A-s\i) = 0, then this specific 
S has no effect on the bound of this theorem as the zero result 
will be discarded by the min* function in |20). 



Case 2: Otherwise, 2~2ies\v P ernl (As\i) > 0. This 
condition ensures that there are no all-zero rows in A$ and 
hence none in Hs(x). However, we know that every J x J 
submatrix of tlg(x) has a zero determinant for the specific S 
that generates c(x) = 0. There are two ways for this to occur. 
One row of Hg(x) is linearly dependent on the other rows or 
two or more rows contain sufficient zeros and zero divisors to 
force the determinant to zero. 

We analyze this case further by setting aside the tth row of 
Hs(x), ht t s[x), preferring a row that is linearly dependent 
on the other rows. If there are no linear dependent rows 
in Hs(x), we choose a row that contains zero divisors. 
We form a new matrix with the remaining J — 1 rows 
of H(x), called H'(x), and another with the corresponding 
J — 1 rows of A, called A'. Because of the assumption that 
2~2ieS\v P erm ( A S\i) > 0' there must be at least one index 
I £ S\V, such that perm (A^wj > 0. The cofactor expansion 



of this term along row t contains a term 

a t,i* •P erm ( A (5V)v) > °> 



(22) 



for some i* £ S\i, where positive value at,i* is the entry in 
the tth row and i*th column of A. Let S* = S\i*. 

To proceed with this case, we assume that H'^. (x) contains 
at least one (J — 1) x (J — 1) submatrix with nonzero 
permanent. If this isn't true, we repeat the row removal process 
above, which may need to be repeated several times. In the 
extreme, H(x) may be reduced to one remaining row. Then 
applying Lemma [T] with H(x) = H'(x) and S = S*, we 
generate a nonzero vector, c*(x), with components 



■(x) 



(Hjf. Xi 0r) 



(23) 



perm | I : ; (.r i j if i £ S* 
otherwise. 

The proof of Lemma [T] implies that H' (x)c* (x) T — T . 
Multiplying the removed row of the parity-check matrix by 
the vector c*(x) yields 

h t (x)c*(x) T = ^2 h t,i[x) ■ c*(x) 
ie[L] 



E 



ht,i(x) perm 



H' s , v (z) 



(24) 



= perm (H s » (x)) = 0, 

since all J x J submatrices of Hs(x) were assumed to have 
a zero permanent. Therefore the nonzero vector c*(x) is a 
codeword in C. 

By puncturing c*(x) we generate another nonzero vector, 
c'*(x), which is a codeword in C. The Hamming weight of 
this codeword satisfies the upper bound 

wt(c'*(a;))= wt (penn(H^. v (a;); 

ies*\v 

^2 P erm ( A 5 



< 



i£S*\P 



< 



perm 

E E at 'J ' perm 
es*\Vjes\i 

perm ( A. 



(S\i)\j 



(25) 



- E 

ies»\v 

< ^2 perm I 

ies\v 

The second line follows from Lemma [2] The third line of ( |25j ) 
has used the fact that (22 1 implies that a t A * > 1. The third line 



( A *V) 



contains a subset of the nonnegative terms of the fourth line. 
The fourth line contains the sum of the cofactor expansions 
of each addend of the fifth line. 

Now, we can state that even a set S which generates the 
all-zero codeword in Lemma [3] provides a valid upper bound 
on the minimum distance of the punctured code 



(C) < J2 p erm ( A ^v) > °- 

i<£S\V 
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V. New Tighter Bounds On Minimum Distance 

Examining the AR4JA protomatrices for rate-2 /3 in ^ and 
rate-4/5 in we see cases where the selection of J + 1 = 4 
columns of the weight matrix A will produce a submatrix As 
containing an all-zero top row. This particular selection of S 
produces the all-zero codeword by the codeword construction 
of Lemmas [T] and [3] and, thus, will have no effect on the upper 
bounds of Theorem [3] and |4] We can improve those bounds 
by finding nonzero codewords after row elimination, as in the 
proof of Theorem |4] 

In the interest of brevity, we will state the following theo- 
rems in a way that applies to both unpunctured and punctured 
codes. In the unpunctured case, it is understood that the set V 
is empty. 

Lemma 5. Let C be a QC code constructed by optionally 
puncturing subblocks of QC code C, defined by the polynomial 
parity-check matrix H(a;) G [F2[ie]/ (x n — l)] . Let the 
subblocks of code C indexed by set V,Vd [L], be punctured, 
while maintaining the dimensionality of the code. Let H (x) be 
a submatrix ofH(x) with rows h t (x), t G T C [J], removed. 
Let S be a subset of [L] of size J + 1 — \T\, such that 



fw (S, t) = perm 



Hfe(x) 



ovi g r. 



(26) 



Let the length-L vector, c'(x) = (c (x) , c[ (x) , . . . ,c' L _ 1 (x)^j, 
with c'^x) G |F 2 [a;]/ (x N — 1^ . ip}, be defined by: 



perm(H^(x)) ifieS\V 
c 'i( x ) = {ip ifieV 
otherwise. 



(27) 



Then c'(x) is a codeword in C. 
Proof: We consider two cases. 

Case 1: If V = (the code is unpunctured), then C = C. 
The product of every retained row of H(.x), denoted hj(x), 
with the vector c'(x) is 



h,(*)c'0r) 3 



E 

ies 



hj,i(x)peTm (n' sv (x) 



= perm 
= det 



h jlS (x) 
U' s (x) 

h jtS (x) 



(28) 



= 0Vj£T, 



since the determinant expression contains a repeated row. 
Next, the product between every removed row of H(x), 
denoted h t (x), and the vector c'{x) is 



h t (x)c'(x) q 



perm 



h t , s (x) 



OV^GT, (29) 



because the permanent was assumed to be zero in ( |26| ). Since 
all rows of the original polynomial parity-check matrix H(x) 
have been accounted for, H(x)c' {x) T — T and c'(x) is a 
codeword in C = C . 



Case 2: If the code C is punctured, let the length- 
L vector, c(x) = (cq(x), Ci(x), . . . , Cl-i_(x)), with Ci(x) G 
¥ 2 [x}/ {x N - 1), be defined by: 



Ci{x) 



perm (H^z 




if ieS 
otherwise. 



(30) 



The proof follows by noting that c'(x) is obtained by punctur- 
ing subblocks indexed by V from the unpunctured codeword 
c(x), above. Since c(x) is in C by Case 1, c'(x) is in C . ■ 
Not only does Lemma [5] remove all-zero rows, it helps 
produce lower weight codewords in more general conditions 
as the following example shows. 

Example 1. Let the following be the columns of the polyno- 
mial parity-check matrix H(x) indexed by set S 



Hs(aO 












x a 


x b 


x c 


x a 


x b 


x d 



(31) 



First, we perform single row removal on Hs(x), since 
perm (Hs(a;)) = as required for |T| = 1. This generates 
all-zero codewords and the codeword subvector 

c s (x) = (x b+d + x b+c , x a+d + x a+c , 0) mod (x N + l) . 

However, looking deeper, Lemma [5] will let us delete two sub- 
rows when the column subset is {0, 1}, producing the obvious 
codeword subvector Cg(x) = (x b , x a , 0), when T — {0, 1} or 
{0,2}. 

Theorem 5. Let C be a QC code constructed by optionally 
puncturing subblocks of QC code C, defined by the polynomial 
parity-check matrix H(x) G [F2[a;]/(x Ar — l)] . Let the 
subblocks of code C indexed by set V, V C [L], be punctured. 
Let H'(a;) be a submatrix ofH(x) with rows h t (x), t G T C 
[J], removed. Let S be a subset of [L] of size J + 1 — |7"|, 
such that fn> (<S, t) = V t G T. Then the minimum Hamming 
distance of C satisfies the upper bound 



drnw. (C ) < mm > wt 

S.T ^— ' 

ies\v 



(perm(H' sv (x))). (32) 



Proof: The proof mirrors the proof of Theorem [3] with 
the weight of the resulting codeword now given by 



wt(c'(x)) = 



E 

ies\v 



wt (perm [H' sv {x) 



(33) 



Note that the minimization in (32i is done jointly over the 
removal of every possible set of rows T and every set of 
retained columns S such that |«S| + \T\ = J + 1 and (26 1 
holds. In the case of single row removal (\T\ = 1), any row 
may be removed, as our requirement ([26]) degenerates to the 



condition perm (Hs(x)) = 0, which is independent of the row 
selected for removal. For multiple row removal, the conditions 
are more complex to evaluate as each row in the set to be 
removed must be tested individually in ( [26] ). 

The exact conditions that allow for row removal in ( |26] i on 
the polynomial parity check matrix for a specific expansion 
factor cannot be duplicated on the nonnegative weight 



submitted to IEEE TRANSACTIONS ON INFORMATION THEORY VERSION: JANUARY 1 1. 2012 



matrix A which is independent of N. Thus, the following 
theorem uses the stricter condition that the sub-row a^s is 
all-zero before removal. 



Theorem 6. Let C be a QC code constructed by optionally 
puncturing subblocks of QC code C, defined by the polynomial 
parity-check matrix H(x) G [F2[x]/ (x n — l)] and let 
A = wt (H(x)). Let the subblocks of code C indexed by set V, 
V C [L], be punctured, while maintaining the dimensionality 
of the code. Let H'(a;) be a submatrix of H(x) with rows 
h t (x), te7"C [J], removed. Let S be a subset of [L] of size 
J + 1 — \T\, such that the sub-rows a t $ = Vi G T, and 
let A! be a submatrix of A with rows a t , t G T, removed. 
Then the minimum Hamming distance of C satisfies the upper 
bound 

~ perm (A' A (34) 



t(C') < min* £ 

' i£S\V 



Proof: The |T| sub-rows of the weight matrix A to be 
removed are all-zero (i.e., & t $ = Vi G T) if and only 
if the corresponding sub-rows of the polynomial parity-check 
matrix H(x) are all-zero (h t s = Vt e T). The latter zero 



condition implies that ( 26 1 holds and we may apply Lemma 
[5] with this S and T and construct a codeword in code C. By 
Theorem |5j the weight of the resulting codeword c'(x) is 



wt(c'(x)) 



wt 



< 



E 

ies\v 

perm 

ies\v 



(perm (H^ v (a:) 



Ai 



s\, 



(35) 



where Lemma [2] is applied in the second line. Once again, we 
must justify the use of the min* operator in regards to the 
all-zero codewords discarded in Theorem [5] We will take a 
specific set S such that every (J — |T|) x (J— |T|) submatrix 
of H' s (x) has a zero permanent. 

Case 1: If 2~2ies\v P erm ( A !s\i) = ^> then this specific 
S has no effect on the bound of this theorem as the zero result 



is discarded by the min function in (34i. 

Case 2: Otherwise, 2~2ies\v P erm ( A s\i) > 0- We set 
aside the jth row of Hs(x), hj t s(x) with j f T, preferring a 
row that is linearly dependent on the other rows. If there are 
no linear dependent rows in H' s (x), we choose a row that con- 
tains zero divisors. We form a new matrix with the remaining 
J — 1 — \T\ rows of H'(x), called H"(x), and another with the 
corresponding J— 1 — |7~| rows of A', called A". Because of 
the assumption that 2~2ies\v P erm (-^-5\i) > 0, there must be 
at least one index i € S\V such that pcrm(A^y-) > 0. The 
cofactor expansion of this positive term along row j contains 
the term 

( A (W)\i*) > 0> (36) 



perm 1 



for some i* € S\i, where positive a^,%* is the entry in the jth 
row and z*th column of A. Let S* = S\i*. 

To proceed with this case, we assume that (x) contains 
at least one (J— 1 — \T\) x (J— 1 — \T\) submatrix with nonzero 
permanent. If this isn't true, we repeat the row removal process 
above. Then applying Lemma[T] with H(x) = H"(x) and S = 



S*, we generate a nonzero vector, c*(x), with components 
'perm (w^ixfj if i G S* 



4 (x) 



otherwise. 



(37) 



The proof of Lemma[T]implies that H" (x)c* (x) T = T . Next, 
the product of the jth row hj(x) that was set aside and c*(x) 
is 

h A x ) c *( x ) T = J2 h jA x )p eTm (H'i.v^x)) 

ies* V ' (38) 

= perm(H^»(a;)) = 0. 

Equation ( |3~8] l follows by noting that the first line is the cofac- 
tor expansion of the second line and the permanent of H's« (x) 
is zero since it was assumed that all (J — |T|) x ( J — |T|) 
submatrices of H' s (x) have a zero permanent. For the rows 
originally removed, we have 



h t (x)c*(x) T =J2°- P erm (Hs-vW) = °> 



(39) 



ies* 



for all t G T, since we noted Yi t ,s{ x ) = at the start. We 
have now multiplied the vector c*(x) by every row of the 
polynomial parity-check H(x). Therefore, H(a;)c* (x) T = T 
and the nonzero vector c* (x) is a codeword in C. 

By optionally puncturing c*(x) we generate another 
nonzero vector, c'*(x), which is a codeword in C. The 
Hamming distance of this codeword satisfies the upper bound 



wt(c'*(a;)) < 



perm 



(aS. V( ) 



<a<j,i*- E P erm ( A s-vJ (40) 

ieS'\v 

< P Crm ( A 5v) 5 

ies\v 

where the justification parallels that of Q25) . Finally, from ( |40| ) 
we can see that this particular S, which produces an all-zero 
codeword in Lemma [5] still produces a valid upper bound on 
the minimum distance of the code, 



perm 



> 0. 



(41) 



Example 2. The sample weight matrix A shown below will 
demonstrate the benefit of Theorem |6l 



A = 



(42) 



Treating the code as unpunctured, Theorem [4] produces a 
minimum distance upper bound of +00, since all 3 x 3 
submatrices of A have a zero permanent. Theorem [6] produces 
a much tighter bound of 3, when A^ = [l 2] . 

Example 3. The following weight matrix A is selected due 
to its similarities to the AR4JA rate- 1/2 protomatrix in Q: 

"0 3 3" 

110 13. (43) 
1 2 2 1 
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TABLE I 

Minimum Distance of AR4JA Protomatrices Stage 2 



TABLE II 

DISTANCE OF CCSDS AR4JA PARITY-CHECK MATRIX 



Code 


Upper Bounds by 


Num. of S sets of 


Code 


Minimum Distance 


Stopping Distance 


Rate r 


Theorems ^ and Q 


size J + 1 in [L] 


Rate 


U.B. by Searching 


U.B. by Searching 


1/2 


66 


7.8 x 10 4 


r 


k = 1024 


k = 4096 


k = 1024 


k = 4096 


2/3 


58 


3.7 x 10 7 


1/2 


52 


63 


50 


63 a 


4/5 


56 a 


5.2 x 10 10 


2/3 


26 




23 


62 b 


"Computations are not exhaustive in S sets due to complexity. 


4/5 


13 


27 


11 


25 



Treating the code as unpunctured, Theorem [4] produces a 
minimum distance upper bound of 30, while Theorem [6] 
produces a substantially tighter upper bound of 10. The reason 
is that Theorem [4] produces distance bounds only with the 
relatively strong contributions of the 3's in the top row of A. 
Theorem [6] will remove the top row in one of its formations 
of which yields a tighter bound. 

VI. QC Expansion of AR4JA 

A direct QC expansion of the AR4JA protographs shown 
in Figs. |4] and [5] will generate a QC-LDPC code. Applying 
Theorems [4] and [6] to the AR4JA protomatrices of Q-Q 
yields upper bounds on the minimum distance of 10 for 
all code rates, independent of block length. As a Hamming 
distance of 10 is rather small for the long block lengths 
desired, a more involved expansion process is of interest. 

The AR4JA codes defined in the experimental CCSDS 
standard J7) use a two-step expansion process. After a first 
cyclic expansion ("lifting") by a factor of 4, a new larger type- 
I weight matrix is obtained as shown in ( 44 1 for rate 1/2. A 



type-I weight matrix is one that contains only ones and zeros 
— meaning that the associated graph does not have parallel 
edges. 



A = 



-o 























1 























1 








1 - 





























1 

















{) 
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1 
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1 
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1 
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1 
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1 























1 








1 
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1 








1 











1 























1 
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1 
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1 
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1 


1 


1 





1 

















1 


1 














1 


1 





{) 


1 





{) 








1 








1 








1 

















1 


1 








1 














1 





1 


1 


























1 


1 








1 





-0 








1 





1 


1 

















1 








1 











1 - 



(44) 



According to the CCSDS standard, the weight matrices so 
obtained, such as (44) , are expanded in a second step cyclic 
expansion to create QC-LDPC codes with three block lengths, 
corresponding to k = 1024, 4096, and 16384 information 
bits. For example, cyclically expanding ( |44) by a factor of 
N = 128 yields the (n = 2048, k = 1024) AR4JA code after 
puncturing. In this final expansion, the binary parity-check 



matrix H is constructed by replacing each 1 entry of (44 1 by a 
cyclic permutation matrix selected by a variation on the ACE 
algorithm. These codes are QC with a subblock size equal 
to the second step expansion factor (e.g., N — 128). In other 
words, the two-step process is not equivalent to any single step 
cyclic expansion. With this in mind, the new protomatrices 



such as (44 1 should be used to compute the QC distance 
bounds described here for proper application to the CCSDS 
AR4JA codes. Those results are shown in Table [I] 



"The smallest stopping set found was a codeword. 
b Beyond the upper bound shown in Tablejlj 



TABLE III 

Weight Spectrum for CCSDS AR4JA Rate-4/5, k = 1024 using 
Search Parameters: I = 150 and T = 10 





Hamming Weight 


Search 
Time 


13 


14 


15 


16 


Num. of 
Codewords 


32 


256 


128 


400 


1.5 hrs 



TABLE IV 

Weight Spectrum for CCSDS AR4JA Rate-4/5, k = 4096 using 
Search Parameters: I = 300 and T = 19 





Hamming Weight 


Search 




27 


28 


29 


Time 


Num. of 










Codewords 


128 








82 hrs 



TABLE V 

Weight Spectrum for CCSDS AR4JA Rate-2/3, k = 1024 using 
Search Parameters: I = 300 and T = 19 





Hamming. 


Weig 


ht 


Search 




26 


29 


31 


32 


Time 


Num. of 












Codewords 


64 


128 


64 


64 


5 hrs 



We have worked to reduce the computation time to produce 
the results of Table [I] For larger weight matrices, if we assume 
that calculations are dominated by the computation time tj 
for the J x J permanent, then the total time to evaluate 
Theorem 4 is tj(J + The final term ( J+1 ) is the 

number or S sets in weight matrix A and is shown in the 
right column of Table [I] for AR4JA. For the computations of 
interest, with J = 12, we built a simple recursive routine with 
ti2 = 44/zs for computing sparse permanents (as measured 
on a 2.6 GHz CPU). Thus, we estimate that computations for 
rate-1/2 in Table [i] will take 44s while we measure a run-time 
of 53s, which includes bookkeeping and set manipulations. 
Unfortunately, we estimate that the rate-4/5 result will require 
344 days to completely evaluate Theorem [4] so we selectively 
direct the computations towards the weaker components to 
yield the result shown in Table [I] Furthermore, our additional 
efforts to recompute the minimum distance bounds using the 
row elimination logic of Theorem|6]did not yield tighter results 
with the AR4JA weight matrices. 
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Fig. 6. Minimum distance bounds vs. block length for rate-1/2 AR4JA. 



VII. Distance Spectrum Found by Search 



Several papers, including | |15] , fl6) , have described search 
techniques to find the minimum distance and/or stopping sets 
of general LDPC codes. To quickly confirm our bounds on 
minimum Hamming distance, we utilize the error impulse and 
decoding algorithm of fl3] for a non-exhaustive search for 
small stopping sets, noting which stopping sets are codewords. 
We modify the algorithm to take advantage of QC symmetry 
by skipping impulse combinations which are identical after 
cyclically shifting every subblock. We also extend the pa- 
rameters (iterations I and maximum threshold T) and erase 
the punctured symbols before the erasure decoding step, in 
addition to the symbols already erased by the algorithm. Our 
upper bounds on minimum distance and stopping distance by 
searching are summarized in Table [IT] 

The selected algorithm is demonstrated on rate-1/2 codes up 
to a minimum distance of 19 in [ [15] , and the weight spectrum 
results are listed only out to a Hamming weight of 25. Thus, 
we may be pushing this algorithm beyond its effective range 
in several cases. Further searching may turn up lower weight 
codewords and stopping sets. 

When a codeword of a QC-LDPC code is located by the 
search technique, it is an indication of a group of codewords 
with the same properties due to the QC nature of the code. For 
instance, the low weight codewords of the rate-4/5, k — 1024 
CCSDS AR4JA code generally occur as a set of N = 32 
cyclically-shifted versions of a base codeword. On occasion, 
the cyclically-shifted codeword returns to the base codeword 
after a fraction of N shifts when all subblocks are periodic 
with a common period. Accounting for this, all unique low- 
weight codewords found via the search have been tabulated 
for several of the AR4JA parity-check matrices given in 
the CCSDS standard in Tables iHlffVl As search time grows 
substantially with increased block length and success rate 
drops with the increased minimum distance achieved at higher 
block sizes, the results are incomplete. 



VIII. Summary of AR4JA 

Fig. [6] summarizes for rate-1/2 AR4IA codes our minimum 
distance results as a function of the block length n. First, 
our QC upper bound of 66 from Theorems [4] and [6] based on 
the weight matrix is plotted across the figure. Indicated by 
'X', the minimum distance upper bounds of Theorem [3] based 
on the QC polynomial parity-check matrices are found to be 
62, 63, and 65 for the block lengths 2048, 8192, and 32768, 
respectively. Finally, from Table |IIJ the weights of codewords 
found by searching are designated by the symbol 'O.' Note 
that the specific codewords found are consistent with the more 
general upper bound result of 66 for all block lengths. 

Divsalar et al. showed that most codes in the ensemble 
of protograph-based codes with certain properties, which the 
AR4JA protograph possesses, have minimum distance linearly 
increasing with block length (5J, (6). Specifically, by upper 
bounding the ensemble average weight enumerator, they were 
able to prove that Pr {<i m i n < <5 m i n n} — > exponentially fast 
as n —> oo. The value <5 m j n = 0.015 has been computed for 
the rate-1/2 AR4JA-based ensemble of codes and used to plot 
the linear minimum distance relationship of Fig. [6] Our fixed 
upper bound on minimum distance applies only to protograph 
expansions that use binary circulant matrices and does not rest 
on probabilistic arguments. Thus, for the QC-AR4JA codes 
that appear in the standard, our bounds dominate after the 
lines cross at n = 4400 bits. 

If we examine the probability that a random expansion is 
quasi-cyclic, we shed some light on the situation. Consider 



the expansion of each 1 entry in a protomatrix such as (44i, 
by a factor N. There are N cyclic permutation matrices to 
choose from and AH general permutation matrices. Thus, a 
randomly chosen permutation matrix has only a probability of 
1 / ( N— 1) ! of being cyclic. This probability goes to zero super- 
exponentially fast. Since the cyclic class of expansions is such 
a small fraction of the ensemble of all possible expansions, 
one cannot claim with certainty that Divsalar's probabilistic 
bounds apply to the cyclic class. 

IX. Girth of AR4JA 
The girth of a code denotes the length of the shortest cycle 



in its Tanner graph. Table VI summarizes our findings on the 
girth of the AR4JA codes. The girths of the standardized codes 
(7) are shown in the table without parentheses for each block 
length and code rate. Next, from the protographs themselves, 
the neighborhood of any node can be diagrammed as a tree. 
We can measure how tall this tree is at a given number of 
nodes corresponding to the specified block length |T], fl4) . 
We do this for each node type in the protograph and select 
the smallest as an upper bound on girth that would apply to any 
possible expansion method. These upper bounds are denoted 



in parentheses in Table VI 



We also apply upper bounds to the girth based on quasi- 
cyclic expansion properties. Since the AR4JA protomatrices 
of Q-Q all contain the element 3, their derived graph by 
QC expansion is limited to girth < 6 flO) . However, since 
the codes in the CCSDS standard use a two-step expansion, 



we must examine the intermediate protomatrices such as (44 > 
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TABLE VI 
Girth of the CCSDS AR4JA codes 



References 



Information 


Measured Girth (Upper Bound) 


Bits k 


r = 1/2 


r = 2/3 


r = 4/5 


1,024 


6 (12) 


4 (10) 


4(8) 


4,096 


8 (14) 


6 (10) 


4 (10) 


16, 384 


10 (16) 


6 (12) 


4 (10) 


QC Limit of 








Protomatrix 


12 


12 


12 



(#) Denotes upper bound computed by tree method Jlj, fl4j 

We find that they contain the submatrix \\ \ \ ] T at every code 
rate. This limits the QC expansion to girth < 12, independent 
of block length, as shown in the final row of Table VI [17|. 

This section has shown that the two-step expansion ap- 
proach was essential for achieving girths beyond 6 for larger 
block-lengths of rate-1/2 AR4JA. 

X. Conclusion 

This work has extended the distance bounds of flO) to 
punctured QC-LDPC codes. We have also tightened those 
distance bounds in cases that are relevant to protomatrices 
that contain many zeros. 

Next we evaluated the minimum distance upper bounds 
for the AR4JA codes specified in the CCSDS experimental 
standard for deep space. We have shown that the two-step 
expansion approach was critical for achieving reasonably high 
minimum distance for these codes in QC-LDPC form. We have 
shown that the minimum Hamming distance of the standard- 
ized QC AR4JA codes does not grow with block length, even 
though the asymptotic ensemble minimum distance of AR4JA 
codes does grow linearly in block length (3), (6). Nevertheless, 
the minimum distance of the CCSDS codes is likely high 
enough for practical purposes. 

The bounds developed here and in [ 10 1 can be useful tools 
in evaluating future QC-LDPC code designs, both punctured 
and unpunctured. 
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